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Abstract 



Helium atoms are strongly attracted to the interstitial channels within a bun- 
dle of carbon nanotubes. The strong corrugation of the axial potential within 
a channel can produce a lattice gas system where the weak mutual attraction 
between atoms in neighboring channels of a bundle induces condensation into 
a remarkably anisotropic phase with very low binding energy. We estimate 
the binding energy and critical temperature for ^He in this novel quasi-one- 
dimensional condensed state. At low temperatures, the specific heat of the 
adsorbate phase (fewer than 2% of the total number of atoms) greatly exceeds 
that of the host material. 



Low temperature research on Helium was 
initially stimulated by the challenge of deter- 
mining the condensation temperature of bulk 
He [Q. In recent decades, two-dimensional 
He films, in which the superfiuid transition 
differs qualitatively from that of the bulk 
0, have been particularly intriguing. While 
once of only academic interest Helium 
in one-dimensional or quasi-one- dimensional 
systems has received increased attention re- 
cently since the realization that such systems 
can be created in the laboratory. He atoms 
are very strongly bound within the hexago- 
nal lattice of narrow interstitial channels be- 
tween tubes within the triangular lattice of a 
bundle of carbon nanotubes [|7HlO[]- Within 
this very narrow channel, the transverse de- 
grees of freedom are frozen out even at rela- 
tively high temperatures of ~50 K. The bind- 
ing energy per atom, 340 K, is the highest 
known for He, almost twice that calculated 
for He within individual nanotubes [|TT| and 



2.4 times higher than that on the basal plane 
of graphite |]T2|,|T3[. It exceeds the ground 
state binding energy of bulk liquid ^He by 



nearly fifty times |14 



Here we describe how the strong axial 
confinement of the He wavefunctions within 
a single channel can produce a direct exper- 
imental realization of a lattice gas model, 
wherein the weak coupling between atoms 
in neighboring channels induces a finite- 
temperature transition into a remarkably 
anisotropic and extremely weakly bound con- 
densed state. First we present a localized 
model wherein the Helium resides in peri- 
odic array of relatively deep potential wells; 
this "bumpy channel" approximation is sup- 
ported by single-particle Helium band struc- 
ture calculations. For comparison, we also 
describe a delocalized model which assumes 
translational invariance within each chan- 
nel (a "smooth channel" approximation). 
The large difference between the models in 
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the transition temperature to the condensed 
phase demonstrates the importance of the 
external potential in controlling the within- 
channel He-He interaction. 

In all calculations we assume that the He- 
He interaction is unaffected by the carbon 
environment, an approximation that omits 
screening by both phonons and electrons. On 
planar graphite the electrodynamic screening 
of the van der Waals interaction is most im- 
portant and reduces the well depth of the He- 
He pair potential by ~ 10% for a monolayer 
film The smaller He-C separation in 

the interstitial channel of a nanotube bun- 
dle should yield a somewhat larger effect; the 
omission of this screening implies a moder- 
ate overestimate in the binding energies de- 
scribed below. 

The external potential that a Helium 
atom experiences due to the Carbon envi- 
ronment of the interstitial channel has a sig- 
nificant corrugation. Using a C-He pair po- 
tential to model this interaction 0, a band 
structure calculation of an isolated ^He atom 
within an interstitial channel of a (18,0) |]T^ 
tube lattice yields a purely one-dimensional 
dispersion with a very large enhancement of 
the lowest-band effective mass m* above the 
bare mass m: m*/m ^ 18, with a band- 
width of 0.18 K (18 times smaller than the 
free-particle bandwidth) |jl8|. For a regular 
lattice of tubes the interchannel tunnelling 
is negligible. Fig. 1 shows contours of 
constant probability density for the ground 
state wavefunction at = in the single- 
particle band structure. The coupling be- 
tween sites is so weak that even at moder- 
ately high temperatures the primary means 
of intersite motion along a channel is single- 
particle atomic tunnelling. Due to the het- 
erogeneity of currently accessible nanotube 
systems (i.e. mixtures of tubes with different 
wrapping angles and diameters), this result 
should be treated qualitatively, as a demon- 
stration that the atomic He states are well- 



localized axially within the interstitial chan- 
nel. Such a small bandwidth implies well- 
confined single-particle wavefunctions which, 
in this particular geometry, occupy a regu- 
lar lattice with a separation of 4.2 A be- 
tween sites. (In the many-body problem, 
the correlation due to hard-core interactions 
with neighbors could induce a further local- 
ization). Note that the single-particle calcu- 
lation described above correctly predicted the 
binding energy per atom, as recently mea- 
sured experimentally |T0[. The bandwidth 
could be controlled by changing the distribu- 
tion of nanotube diameters [W\ or through 
external pressure, possibly inducing a quan- 
tum phase transition. 

For this lattice of localized Helium sites, 
the natural description is a lattice gas model 
wherein the statistical degrees of freedom at 
low energy are the occupation indices of the 
or 1 



I.e. 



sites, 
ping energy 



The intersite hop- 
i.e. ~0.1 K) is significantly 
lower than the potential energy of interac- 
tion between atoms on neighboring sites in a 
channel (~0.5 K) [|2T]]. Multiple occupancy is 
excluded by the hard-core repulsive interac- 
tion between atoms (estimated to impose an 
energetic cost of at least 100 K). Intra-site 
excitations also involve high energy scales, 
which are irrelevant at low T. The calculation 



uses results previously obtained ||2^ , [23[| for 
an anisotropic simple cubic Ising model with 
an interaction strength between neighbor- 
ing spins along the z axis (i.e. within the 
same channel) and a transverse interaction 
Jt = c Jz- For the He-nanotube system, c 
is very small. When c < 0.1, the transi- 
tion temperature is well approximated by the 
asymptotic formula [0,23] 



1) 



J, ln(l/c) - ln[ln(l/c)] 

Here we have Jz = \V{a)\/4:, where V{a) 
is the equilibrium interatomic interaction at 
the intersite separation and the factor of 1/4 
arises from the familiar transformation from 
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the Ising model to the lattice gas. For Helium 
with an intrachannel site spacing of a = 4.2 
A, Jz = 0.5 K. As the present lattice is honey- 
comb rather than simple cubic and as second 
neighbors are not included in Eqn. |I], we esti- 
mate that the transverse interaction strength 
Jt should be renormalized by roughly a fac- 
tor of 3 X 3/4, which includes a factor of ~3 
for the second neighbors in adjacent channels 
and a factor of 3/4 for the reduced coordina- 
tion of the lattice. For Helium with inter- 
channel spacing d = 9.8 A we obtain Jt ^ 7 
mK, and c ~ 0.015. Eqn. ([l|) then yields 



0.7 J. 



(2) 



The transition temperature for condensation 
in this lattice-gas Helium model is then ~ 
0.36 K. Reducing c by even a factor of 5 
would change Tc by less than 25%, illustrat- 
ing the insensitivity of Tc to the transverse 
interaction. Variation in Tc through changes 
in will be minor so long as this lattice 

gas model remains valid. The binding en- 
ergy of this system is lower than that of any 
other stable condensed atomic system, to our 
knowledge. 

Although this system falls in the three di- 
mensional Ising model class, the quasi-one- 
dimensional character of the system mani- 
fests itself in the specific heat above the tran- 
sition. For small c, the specific heat above the 
transition p3| (shown in Fig. 2 at the critical 
site occupancy of 1/2) closely follows the one 
dimensional result, 



C(T)/(kBN) 



knT 



sech 



J. 



kRT 
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(3) 



The one dimensional specific heat has a max- 
imum of ~ 0.44 near T/Jz = 0.83, slightly 
above the transition. We take this tempera- 
ture as a convenient reference point for com- 
parison to other systems. For d = 9.8 A 
and a = 4.2 A, the resulting adsorbed par- 
ticle density implies a specific heat of ~ 4 



mJ/gK (normalized to the mass of the car- 
bon host) at T 0.83J^. This contribution 
to the specific heat should be observable: it 
greatly exceeds the background specific heat 
of the host material. For example, measure- 
ments on a sample of single-walled carbon 
nanotubes yield a specific heat of C ~ 0.2 
mJ/gK at 1 K which is decreasing with de- 

m 



creasing temperature 



The specific heat 



of graphite is at least three orders of magni- 



tude smaller in this temperature range p5 



Note also that experimental measurements 
of the specific heat of single-walled nanotube 
bundles significantly exceed theoretical esti- 
mates for the contributions from the nan- 
otube substrate 



the low-energy degrees 
of freedom of adsorbed gases might account 
for this discrepancy. 

Comparison of these results to calcula- 
tions in the smooth channel approximation, 
wherein the external potential is flat, reveals 
the importance of the confinement to defined 
lattice sites in this strongly anisotropic con- 
densed He. Diffusion Monte Carlo calcula- 
tions for a one dimensional assembly of "^He 
atoms in a fiat potential [^,^ yielded a very 
weakly bound state (~ 2 mK per atom) of 
remarkably low density (~ 0.04 A~^). We in- 
troduce the interchannel interactions through 
a variational wave function which is a prod- 
uct of identical states of density p in every 
interstitial channel. The energy shift A due 
to the interchannel interaction between chan- 
nels separated by distance d is 



p / dx V{r') 



(4) 



where r' = \/x^~+lP. A straightforward nu- 
merical integration of (0), using the Aziz in- 



teratomic potential for Helium i^" 
eluding the three nearest and the six next- 
nearest neighbor channels yields A/p = 
—0.228 K A. Because of the rapid decay of V 
with distance, inclusion of more distant chan- 
nels produces no appreciable change. Fig. 
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3 shows how the interchannel interaction in- 
creases the binding energy (2 16 mK) and 
the equihbrium density (0.035 —>■ 0.080 A~^) 
of the hquid state above that in the single- 
channel picture. 

Classical statistical mechanics implies a 
proportionality between the ground state co- 
hesive energy and the critical temperature of 
a given system. Although this law of cor- 
responding states fails for quantum systems, 
a specific class of systems, (i.e., a definite 
De Boer quantum parameter p8|), typically 
has a strong correlation between the criti- 
cal temperature and the ground state cohe- 
sive energy. For example, ^He in three di- 
mensions has a binding energy per particle 
of 7.17 K and a critical temperature 5.2 K. 
In two dimensions, these values are 0.87 K 
and 0.85 K, respectively p9| , ^ . The smooth 
channel model should then condense at ~ 10 
mK. The variational approximation underes- 
timates Tc, whereas the neglect of fluctua- 
tions in this nearly one dimensional system 
overestimates Tc. 

The delocalized smooth channel approx- 
imation yields a much lower transition tem- 
perature (~ 10 mK) than the localized case 
(~ 0.3 K). For the smooth channel, the ki- 
netic energy maintains a large (~ 15 A) 
distance between He atoms within a chan- 
nel, which reduces their interaction energy 
by about 30 times relative to the bumpy 
channel approximation. In the more realis- 
tic bumpy channel model, the kinetic energy 
arises mainly from the curvature of the site- 
localized single-particle wavefunction and is 
not relevant to the condensation tempera- 
ture. The external axial potential forces a 
much smaller He-He separation and thereby 
produces a much higher transition tempera- 
ture to the condensed phase. 

In experimentally produced bundles of 
nanotubes the distribution of nanotube di- 
ameters is rather sharp, but present evi- 
dence suggests a heterogeneous distribution 



of wrapping angles for the graphene layers 
comprising individual tubes. The consequent 
variations in the binding and axial separa- 



tion of adsorbtion sites ||3T[ can strongly af- 
fect the critical phenomena near the transi- 
tion and the nature of the condensed phase, 
but should cause comparatively minor vari- 
ations in the transition temperature of con- 
densation. As Tc is linearly proportional to 
Jz and only weakly dependent on Jf in this 
regime, the main effect on Tc will arise from 
variations in the intra-channel intersite sep- 
aration. Should the spatial extent of axial 
localization be substantially longer in some 
bundle geometries, then double occupancy 
could become important as well. 

In summary, we demonstrate that He- 
lium within the interstitial channels of carbon 
nanotube bundles can condense into a novel 
strongly anisotropic phase wherein the strong 
axial confinement induced by the external po- 
tential of the surrounding nanotubes greatly 
enhances the density and hence the transi- 
tion temperature to the condensed phase. At 
low temperatures the specific heat of the he- 
lium absorbates, which comprise less than 2% 
atomic fraction, is predicted to exceed greatly 
that of the much stiffer background mate- 
rial. While we focus on '^He, this intrigu- 
ing interstitial-channel condensation could be 
observable for other small atoms, for exam- 
ple, H2 and Ne 0. Preferential adsorp- 
tion to the interstitial sites can be guaranteed 
by choosing small adsorbates which energet- 
ically prefer the tighter coordination of the 
interstitial channel over the interior of a tube 
(or by simply using nanotubes with closed 
ends). For ^He, we expect the ordering be- 
havior to be very similar to that of "^He (ex- 
cept for the possibility of a very low tempera- 
ture (~ fiK) magnetic transition in the Fermi 
species), since atomic exchange is negligibly 
small. This situation is reminiscent of the 
epitaxial ordering transition of ^He and ^He 
on the graphite basal plane. In that case. 
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exchange is much higher than here, yet tran- 
sition temperatures of the isotopes differ by 
less than 1% 0,13]. 
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FIGURES 



FIG. 1. Isosurfaces of constant probability density for the k = state in the lowest band for 
He atoms within the interstitial channel of a lattice of (18,0) nanotubes. The axis of the chanel 
is oriented horizontally. The three surfaces correspond to probability densities of lO'', 10~^, and 
10"'^ A~'^. 99% of the probability is within the middle surface and 99.98% lies within the outer 
surface. The axial latttice constant is 4.2 A. The lowest-band states are well-confined to distinct 
lattice sites. 

FIG. 2. Specific heat as a function of reduced temperature at one-half site occupancy in the 
anisotropic lattice-gas He model with = 0.5 K (thin curve). The specific heat diverges at Tc = 
0.36 K and asymptotes to the one dimensional specific heat (dotted curve) at high temperatures. 
The lower (thick) curve shows a theoretical result for the specific heat of the nanotube substrate 
(x 1000) 0. 

FIG. 3. Energy e per ^He atom in an isolated smooth one-dimensional channel (solid line), com- 
pared to the energy per atom for a hexagonal lattice of smooth channels at the same interchannel 
separation as for a carbon nanotube bundle (dashed line). 
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